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GRNDYE 1970 - A PROGRAM TO ESTIMATE CARDIAC OUTPUT 
OFF-LINE FROM AN INDICATOR DILUTION CURVE 


DECUS Program Library Write-up DEC US NO. 8-524 

INTRODUCTION 

A popular technique of estimating cardiac output in laboratory 
animals and humans has been the indicator dilution method. Although 
its accuracy has been questioned * , the Stewart-Hamilton method 0 
with indocyanine green as the indicator, has been most widely used. 
Basically, a known amount of dye is delivered as a bolus into the pul¬ 
monary artery, as aortic blood is sampled for dye concentration. The 
concentration of the dye in the aortic blood is determined by a densi¬ 
tometer and recorded, against time, by an electronic recorder. Cardiac 
output is then calculated from the equation: 

. Amt. injected 

(U....Q - ~ps- Where Q = cardiac output 

J o c(t) dt c(t) = concentration of dye 

One of the chief problems in this method is the effect of recirc¬ 
ulation. Equation (1) assumes that there is no recirculation, i.e., 
that no dye which has flowed past the sampling port is allowed to re¬ 
circulate past the sampler. Such of course is not the case in the 
intact animal or human, and recirculation does occur. Although the 
effects of recirculation can be minimized by the use of catheter 
densitometers with short time lags^, most laboratories employ densitom¬ 
eters which are not placed in the aortic blood stream and resulting 
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dye concentration curves show the effect of recirculation before con¬ 
centration has dropped to a negligible level. The Stewart-Hamilton 
method fits a negative exponential function to the portion of the 
curve, c(t), from just past its maximum to infinity, in order to 
ignore mathematically the effects of recirculation. 

This paper reports a program written in PAL-8 machine language 
which, given calibration data and measurements from an indicator dilu¬ 
tion curve, computes cardiac output by the Stewart-Hamilton method. 

REQUIREMENTS 
PDP-8 with 4K; High Speed Reader. 


METHODS 


1. Obtaining the Dye Curve : 

A catheter is positioned in the main pulmonary artery and a 
second catheter in a peripheral artery or in the aorta of the subject. 

A densitometer is used to measured arterial concentration of dye and 
is connected to a recorder. The densitometer should have a switch 
which gives a known signal to the recorder, as a check on the ampli¬ 
fier settings and should be proven to be linear in its response. This 
check deflection is measured before each curve (and called Curve Check) 
and before calibration (and called Calibration Check). 

Dye curves are obtained by bolus injection of indocyanine 
green from a special volumetric syringe (0.75 to 1.5 ml) into the 
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pulmonary artery. The densitometer is balanced and zeroed immediately 
before the injection. 

Calibration of the system is done as follows: The amount of 
dye injected at each dye curve (same volume, same batch) is mixed with 
10 ml of the subject's blood. 0.4 ml of this mixture is diluted with 
9.6 ml of the subject's blood. An additional 10 ml of the subject's 
blood is drawn through the cuvette, while the densitometer is zeroed 
and balanced. The diluted dye is then drawn through the cuvette and 
the resultant deflection (in mm) is called Calibration Factor. 

Data are then obtained from each recording: i) from each sub¬ 
ject, the Calibration Check and Calibration Factor (in mm) and ii) from 
each dye curve, the Curve Check and the deflection of the curve (in mm), 
every 0.5 sec. Curve measurements are taken from the time of departure 
of the curve from the baseline to its return near baseline (usually 6-10 
secs, later—see Fig. 1). These data are presented to the program 
for computation. 

In the event that no electronic check on the amplifier setting 
is provided by the densitometer, both Curve Check and Calibration Check 
can be called equal to 1. 

2. Preparation of Data Tape : 

EDITOR (DEC-08-ESAB) is used to prepare the data tape in ASC ++ 
format. For each set of curves the following are given as integers: 

i) Calibration Check (carriage-return-line-feed) 
ii) Calibration Factor (carriage-return-line-feed) 





Then the data for each curve in the set are given seriatim 
in integer format: 

i) Curve Check (carriage-return-line-feed) 

ii) D1; D2; D3; D4DN; (carriage-return-line-feed) 
where the D's are deflections in mm, and semicolons are used as termin¬ 
ators for each deflection. These deflections must be on one line of 
the data tape, and there can be no more than 24 such deflections given. 

(The program ignores all deflections subsequent to the 24th). 

When the set of curves is finished, the Curve Check for the 
next line is coded as zero ("0; carriage-return-line-feed"). 

If more sets are to be done, the process above is repeated, 
giving the Calibration Check, Calibration Factor and Data from each 
curve as above. The number of curves in a set must be less than 4096; 
the number of sets is limited only by the paper-tape available. 

When the data tape is complete, two entries of "0; carriage-return- 
line-feed" are appended. 

3. Running the Program : 

The binary tape, GRNDYE 1970, is loaded after Floating Point 
Interpreter #4 (Digital-8-25-F), by use of the binary loader (DEC-08-LBAA). 
The data tape is placed in the High Speed Reader, and the program is started 
at Starting Address = 0200. 







DESCRIPTION OF PROGRAM 


I. Data : 

The program accepts positive numbers less than 4096 on the high 
speed reader with as the terminating character for each number. It 
stores them in buffer, with the octal number in the first of three 
locations and 0's in the succeeding two locations. The program does this 
in a subroutine, called Decrd, which is a modification of DIG-8-28-U-SYM. 

The terminating character. Carriage-return-line-feed, causes exit from 
this subroutine. Decrd counts entries as it reads them (and calls the 
number, N) and checks to make sure the number of entires does not exceed 
24. If so, it ignores subsequent entries. 

The program interprets the first number in buffer as "Calibration 
Check" on the first pass through Decrd; on the second pass, it inter¬ 
prets the first number in buffer as "Calibration Factor". For subse¬ 
quent input, the program interprets the first number in buffer as "Curve 
Check"; deflection measurements are then read and counted on the next 
pass through Decrd. The program later recycles through the "Curve Check" 
and deflection readings until it encounters a Curve Check = 0. 

In the event that the program encounters a "Curve Check" equal to 
0, it reinitializes; so that different "Calibration Check" and "Calibration 
Factor" can be read in on the same tape. Should it encounter a "Calibration 
Check" equal to 0, it terminates. 

II. Calculation of Area under the Dye Curve : 

I. Recognition of "upswing" and "downswing" : 

The program compares each deflection with the preceding one until 








the difference is negative, calling this deflection the maximum. It then 
compares subsequent deflections until it encounters a deflection more 
than 5 MM below the maximum. The ordinal number of this deflection is 
called Tmax (TMCTR). The data up to Tmax are considered to be in the 
"upswing" and subsequent data from the "downswing" of the dye curve. 

2. Calculation of area under upswing : 

This is done by Simpson's Rule, after the program recognizes 
Tmax as even or odd. The area is called Al. 

3. Calculation of area under downswing : 

The Stewart-Hamilton method calls for the fitting of a negative 
exponential curve from Tmax until recirculation occurs. Recirculation 
is recognized by the fact that all subsequent observed points lie above 
the fitted negative exponential curve. 

The program finds the natural logarithm of each deflection from 
Tmax to N. It then fits a straight line by least squares to Tmax and the 
next two points. All subsequent points are checked against theoretic 
points of this semi logarithmic fit. If all lie above the corresponding 
theoretic points, it recognizes TMCTR + 3 as "recirculation" and calcu¬ 
lates the answers. If one actual deflection is less than the theoretic 
deflection, the program recycles and fits to Tmax and the next three 
points. This process is continued until recirculation is recognized or 
until N-l is reached. In the latter case, the program prints "NO RECIRC" 
and computes the answers. 

Area is calculated by finding the integral from Tmax to infinity 
of the negative exponential fit, and is called A2. 





III. Calculation of Cardiac Output ; 

The following formula is used to calculate cardiac output: 

Q = 15*(calib factor)•(curve check/calib check)/(Al + A2) 

IV. Output : 

The program prints calibration check and calibration factor and 
then for each curve types out the following: curve number, curve check, 
and all deflections given to program (up to maximum of 24). It then 
reports Tmax (in seconds), recirculation time (in seconds), the sum 
square (variance) of the semi log least square fit, and Q, the cardiac 
output (in L/min). 

Two output programs are used: a subroutine, called Decprt 
(modified DIG-8-22-U-SYM) and the output controller of floating point 
interpreter #4. Both use the (wait for flag; JMP.-l) routine, as does 
Decrd, and must be modified if the interrupt is turned on. 


V. 



Modifications Easily Accomplished : 

1. To increase maximum number of deflections allowed: 

On source tape, change BUFFR = 717 to BUFFR = 1744; change 
MAXNUM = -31 to MAXNUM = -1211; and assemble source tape. 

2. For paper tape system (program is presently geared to run off 
a disc or dectape system) toggle in a HLT (7402) in location 
223, with the program in core. 

3. To assemble the source tape with a PAL-III or PAL-D assembler, 
remove all "FIXMRI" pseudo-operations from source tape. 
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CHECK OF PROGRAM'S ACCURACY 


Fifty-six indocyanine dye-dilution curves were obtained on 
14 miniature pigs by the methods outlined above. For each curve, the 
area under the "non-recirculation" curve (A1 + A2) was obtained by 
hand-calculation by the Stewart-Hamilton method and from the program, 
GRNDYE. These areas are given in Table 1. 

Results were then compared for agreement. If hand-calibrated 
data are considered as the independent variable and program-calculated 
areas as the dependent variable, then the least squares fitted straight- 
line^ has the form: 

y = 0.999x - 1.880 
with r = 0.9996 (PC.01) 

Figure 2 shows the points of these experiments along with the 
line of identity and demonstrates graphically the excellent agreement 
between computed and hand-calculated results. 

CONCLUSIONS 

An off-line program is presented which calculates cardiac output 
on a PDP-8 computer. This program is written in PAL machine language 
and is usable in a paper tape or disc-monitor system. This program is 
sufficiently accurate, compared with hand calculations, to be useful to 
the laboratory measuring cardiac output by indicator dilution methods. 

It suffers from the limitations of the Stewart-Hamilton method itself. 
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Table I 


HAND CALC. VS. GRNDYE - AREAS 

Hand Hand Hand 

Calc. Comp. Calc. Comp. Calc. Comp. 


436.3 

431.8 

358.4 

367.8 

299.5 

297.1 


442.6 ’ 

443.7 

380.9 

362.3 

370.1 

355.0 


226.0 

226.5 

383.3 

365.1 

339.4 

326.9 


178.6 

178.1 

243.9 

233.6 

194.6 

188.0 


185.9 

190.0 

242.0 

239.4 

231.1 

243.3 

ji 

239.4 

246.4 

316.2 

314.1 

192.0 

178.6 

V 

126.6 

126.7 

439.4 

437.3 

240.7 

239.6 

i 

409.3 

447.0 

291.8 

289.1 

290.2 

288.0 

I 

f 

1 

434.8 

472.0 

398.5 

391.0 


• 

\ 

} 

471.9 

465.5 

217.2 

215.1 



t 

i 

\ 

1 

314.9 

315.2 

210.2 

210.5 

ji * 


i 

332.2 

319.1 

230.7 

227.3 

1 

i 


i! 

* 5 

330.5 

334.5 

i 450.9 

i 

448.4 

i 


I \ 

i! 

248.5 

246.5 

1 

1 437.3 

< 

441.2 



! ! 

212.7 

211.1 

1 218.8 

1 

214.6 



1 , 

219.7 

206.3 

| 419.7 

419.2 

« 

i 


\\\ 
i i 

: j 

470.0 

. 463.5 

| 402.6 

1 

398.9 

« 

t 

) 

1 

J 


t. V 

•*> 

1 

482.7 

472.7 

376.1 

374.1 



r 

399.7 

394.6 

379.3 

378.1 

• 


!jl 

353.0 

352.2 

360.8 

362.9 



>* 

it 

375.0 

368.7 

516.4 

522.1 




491.5 

487.0 

508.9 

% 

505.0 


i 

|l 

i; 

485.0 

471.0 

! 

540.3 

535.0 


j 

i 

• 1 

| 

451.1 

451.9 

1 

| 534.9 

i 

t 

522.3 



it 



1 

I 

10 



i 

l 
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Fig. 1 Typical Dye Curve 







Fig. 2 Hand calculated areas versus areas calculated 
by program. The straight line superimposed 
on the data is the line of identity. 
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APPENDIX I: TYPICAL DATA TAPE 


1201 


1231 
120 ! 

9122147170!79173160!45135127122119!181 


118 ! 

9l27!51!63! 57146! 34125118115113! 121 111 10191 


0! 




APPENDIX II: TYPICAL OUTPUT 


GRNDYE 1970 

CAL CH* 120 CAL P* 123 


CRV NOt 

1 

CM* 

120 

4 


9 

22 

47 


70 

79 

73 

60 

45 


35 

27 

22 

19 

18 




TM* 3 < 

40 SEC 

TR* 5.0 

SEC 



SUM 

SO* 0.1704196E-02 




Q* 

6.573 

L/MIN 





CRV NOt 

9 

2 

27 

CM* 

51 

118 

63 

57 

46 

34 

25 


18 

IS 

13 

12 

11 


10 

9 

TM* 2.5 

SEC 

TR* 4.5 

SEC 




SUM SQ* 0.4476360E-02 


Q= 9.441 L/MIN 


GRNDYE 1970 




